Non- conservative forces and effective temperatures in active polymers 
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We follow the dynamics of an ensemble of interacting self-propelled semi-flexible polymers in 
contact with a thermal bath. We characterize structure and dynamics of the passive system and 
as a function of the motor activity. We find that the fluctuation-dissipation relation allows for the 
definition of an effective temperature that is compatible with the results obtained by using a tracer 
particle as a thermometer. The effective temperature takes a higher value than the temperature 
of the bath when the effect of the motors is not correlated with the structural rearrangements 
they induce. Our data are compatible with a dependence upon the square of the motor strength 
(normalized by the average internal force) and they suggest an intriguing linear dependence on the 
tracer difi'usion constant times the density of the embedding matrix. We show how to use this 
concept to rationalize experimental results and suggest possible innovative research directions. 
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I. INTRODUCTION 

The constituents of active matter, be them particles, 
lines or other, absorb energy from their environment or 
internal fuel tanks and use it to carry out motion. In 
this new type of soft condensed matter energy is par- 
tially transformed into mechanical work and partially dis- 
sipated in the form of heat [l| . The units interact directly 
or through disturbances propagated in the medium. 

In systems of biological interest, conservative forces 
(and thermal fluctuations) are complemented by non- 
conservative forces. Realizations of active matter in biol- 
ogy are thus manifold and exist at different scales. Some 
of them are: bacterial suspensions [2-5], the cytoskele- 
ton in living cells [6-10], or even swarms of different an- 
imals (TTj . 

Clearly enough, active matter is far from equilibrium 
and typically kept in a non-equilibrium steady state. The 
difference between active matter and other driven sys- 
tems, such as sheared fluids, vibrated granular matter 
and driven vortex lattices is that the energy input is lo- 
cated on internal units (e.g. motors) and therefore homo- 
geneously distributed in the sample. In the other driven 
systems just mentioned, the energy input occurs on the 
boundaries of the sample. Moreover, the effect of the mo- 
tors can be dictated by the state of the particle and/or 
its immediate neighborhood and it is not necessarily fixed 
by an external field. 

The dynamics of active matter presents a number of 
interesting features that arc worth mentioning here. Ac- 
tive matter displays out of equilibrium phase transitions 
that may be absent in their passive counterparts. The dy- 
namic states display large scale spatio-temporal dynam- 
ical patterns and depend upon the energy flux and the 
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interactions between their constituents. Active matter 
often exhibits unusual mechanical properties, very large 
responses to small perturbations, and very large fluctua- 
tions - not consistent with the central limit theorem. 

Much theoretical effort has been recently devoted to 
the description of different aspects of these systems. The 
investigation of the self-organization of living microor- 
ganisms is summarized in [12l|. The identification and 
analysis of states with spatial structure, such as bundles, 
vortexes and asters was developed in, e.g., (l3l - [T5| . The 
study of the rheological properties of active particle sus- 
pensions was carried out in [16l4l9{ with the aim of grasp- 
ing which are the mechanical consequences of biological 
activity. 

Surprisingly enough, molecular dynamics computer 
simulations, that have proven to be so helpful to probe 
and understand the dynamics of complex passive sys- 
tems, have not been much employed in the field of active 
matter yet. (For a recent review of the most used compu- 
tational approaches to cell mechanics see Ref. ^].) Stud- 
ies presented so far have been mainly analyticaljbased on 
refined calculations on very stylized models [3, 13, [TTl - [l8l | , 
computer simulations of continuum models, e.g., Lattice 
Boltzmann equations (2l| . and numerical studies of rel- 
atively simple lattice models as the ones pioneered by 
Vicsek and collaborators 0, i, [H [l!] . 

Quite independently, passive systems with complex out 
of equilibrium dynamics have been the focus of much 
attention too. These are systems with competing in- 
teractions between their constituents - and no external 
energy input - that in some conditions - such as suffi- 
ciently low temperature or high density - cannot equili- 
brate with their environment. A regime (or phase) with 
glassy properties then occurs. In the last two decades or 
so, a rather complete mean-field theory that describes, 
up to a certain extent, the super-cooled liquid and deep 
glassy regime, developed ,24-28]. Keywords are "mode- 
coupling theory " and the "random first order phase tran- 
sition scenario" . This has been complemented with ex- 
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tensive molecular dynamic simulations of microscopic 
models with interactions between their constituents that 
are as realistic as possible [13, [13 • Simulations have 
shown agreement with certain aspects of the theoretical 
approach, especially the ones that will be of interest in 
this paper. 

An important outcome of such theoretical modeling 
of glassy systems, later confirmed with numerical sim- 
ulations, is the generation of an effective temperature, 
Teff, in super-cooled liquids driven out of equilibrium 
and glasses in relaxation j2^, [s^l- The concept of ef- 
fective temperatures in out of equilibrium systems, as 
arising from deviations from the fluctuation-dissipation 
theorem (FDT) (sij , appeared in studies of dynamical 
systems [S^l and it was later extended and developed in 
the field of glassy systems [2^ . It was here fully appre- 
ciated that Teff should be understood as a time-scales 
dependent parameter and that it takes a thermodynamic 
meaning in systems with slow dynamics only [2^, [s^ . It 
was later proven that in systems with sufficiently slow 
dynamics - reaching a regime of small entropy produc- 
tion - the same values of the effective temperature are 
obtained with a variety of other measurements. Experi- 
mental measurements of T^ff in different glassy systems, 
including glycerol (34] . s pin-gla sses ^35j , poly mer sam- 
ples [ll] , granular matter |37l 138} , laponite |39j | and other 
colloidal suspensions [i^ have been performed. Very re- 
cently, a study of T^ff in an active colloidal suspension 
of Janus particles under a gravity field appeared in (4]| 
and we shall comment on this work later in the context 
of our results. On the numerical side, the work in is 
especially relevant to our study and we shall also com- 
ment on it below. The validity of FDT in the Vicsek 
model was investigated in [23 ] . In the latter paper a 
dynamic measurement, similar to one of the calculations 
we present below, was performed. 

In [43j we took a first step in the direction of ana- 
lyzing the structure and dynamics of active matter with 
molecular dynamics numerical methods. We studied the 
dynamics of an ensemble of self-propelled particles in 
contact with an equilibrated thermal bath, a reasonable 
model for real but simple active matter, such as bacterial 
colonies. In particular, we analyzed the generation of an 
effective temperature and some of its properties. In this 
article we go beyond this study to get closer to more com- 
plex cases. We introduce a model for biological structures 
made of filamentous semi-flexible polymers 4J] . We ana- 
lyze its structure and dynamics and we study the effective 
temperature by using a variety of independent measure- 
ments. In particular, we show how to measure the effec- 
tive temperature by using techniques that have already 
been exploited in the study of the mechanical properties 
of real biological systems. The idea is to mimic real ex- 
periments which probe mechanical properties of the em- 
bedding matrix by following the dynamics of tracer par- 
ticles (see, among many others, Refs. (45l - l47| V both free 
and driven by external fields, e.g., optical tweezers. We 
demonstrate that these techniques can be useful to obtain 



a direct characterization of the out-of-equilibrium state, 
and allow for the determination of an effective tempera- 
ture. A short summary of some of the results included 
in the present work can be found in Ref. (isj . 

The paper is organized as follows. In Section |ll] we 
give details of the model and the molecular dynamics 
computer simulations performed. Section IIIII presents a 
study of the molecular model in its passive limit in which 
we determine the optimal parameters, in particular the 
polymer length, to be used in the simulations of the ac- 
tive case. In Section |TV] we describe our results on the 
structure and dynamics of the molecular active matter 
model. Section |V] contains a detailed analysis of the ef- 
fective temperature measured by different methods and 
a comparison with other studies in the literature. Our 
entire set of data is discussed and put in more perspec- 
tive in Section IVIl where we also give some conclusions 
and a few hints for future work. 



II. MODEL AND METHODS 

Many complex biological systems are constituted of 
elongated fibrillar structures (filaments) J4]. Although 
one could expect to use the well-established polymer the- 
ory to describe such systems, this theory focuses on the 
limits of completely fiexible or totally rigid polymers |4^, 
[50| and most biologically relevant polymers or filaments 
are in neither of these classes. They rather belong in 
a third intermediate category: the one of semi-fiexible 



filaments with persistence and contour lengths 5l[ of 
comparable magnitude. The contour length of these fila- 
ments is too small for them to form loops and knots, yet 
they are sufficiently flexible to have signiflcant thermal 
bending. These problems are much harder to deal with 
analytically [52j and, surprisingly enough, relatively lit- 
tle numerical work has been done so far. In this paper 
we start fllling this gap by presenting molecular dynam- 
ics simulations of a simple but sufficiently realistic model 
mimicking the behavior of semi-flexible (or semi-rigid) 
polymers in interaction. 



A. Semi-flexible polymer model 

We consider a model of a (passive) semi-flexible poly- 
mer in which the components have the shape of fllaments. 
Each fllament is a linear chain of beads. The optimal 
values of the parameters have been fixed by preliminary 
calculations, detailed in [s^. We briefly motivate them 
here. 

Polymer chains are coarse-grained and each segment 
is formed by a bead. Each bead is in contact with a 
thermal environment that is described by a random noise 
and a viscous drag. The ith monomer pertaining to chain 
labeled a has coordinates r^a and velocity Via. The latter 
evolve following the Langevin equation: 
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TjjQ is a Gaussian white noise representing thermal ag- 
itation with zero mean and variance {f]ia{i)v'jbi^')) ~ 
2^miaT6{t — t')S'^"^SijSab, where fi,^ — 1, . . . ,d and d is 
the dimension of space (c? = 3 in what fohows) . The tem- 
perature of the bath is T and we set /cb = 1 henceforth. 
The term —^rrna'Via is the frictional force. We take into 
account the over-damped character of the dynamics by 
considering a large value of the friction coefficient. The 
numerical integration of Eq. ([1]) generates a time evolu- 
tion that is consistent with the canonic ensemble (NVT) 
in thermodynamic equilibrium. Note that we completely 
disregard the effect of hydrodynamic interactions due to 
internal and/or intermolecular degrees of freedom. We 
consider identical beads, with masses rriia = m and di- 
ameter Via = tq for all i and a. 

The deterministic mechanical conservative force on 
monomer i pertaining to polymer a due to all 
other components is ff„ E^i E^i fyab = 




i.e., the usual Lennard- Jones interaction potential cut 
and shifted to zero at Tc by adding the constant value 
—Uivc). In addition, we consider = correspond- 
ing to the position of the minimum of the pure Lennard- 
Jones potential, and obtain a purely repulsive potential 
(C/(rc) — —e in this case). The Heaviside function 6 
assures that the potential is zero for r > Tc- This par- 
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Connectivity of the chain is assured by harmonic springs 
of elastic constant k and equilibrium length ro, acting 
between nearest neighbor monomers. Chain rigidity is 
controlled by a repulsive Lennard Jones potential be- 
tween next nearest neighbors, with Tc = 2^^^S and again 
U{rc) = Finally, monomers separated by more than 
two bonds interact via Eq. ([3]). The actual choice for 
the parameters fc, ro and S allows one to control the 
flexibility of the chains, ranging from semi-rigid to flex- 
ible chains (ssj . In particular, it is evident that chain 
rigidity can be induced by a mismatch of the zero point 
energy position of the two terms in Eq. (|4]) taken sep- 
arately. Therefore, the relevant parameter is the ratio 



- Sb=i J2j=i ^iaUtotirfj'), where rf^ is the distance be- 
tween monomer i in polymer a and monomer j in poly- 
mer b. N,n and Np are the total number of monomers 
per chain and the total number of polymers respectively. 
All polymers have the same length iV^. We consider the 
interaction potential 

Utot{r) ^ Uintra{r) + Uinter{r), (2) 

in which the two terms are the intramolecular (a =/= b) 
and intermolecular (o = b) contributions, respectively. 
The exact form of the two contributions must be chosen 
in order to expand the region of stability of the liquid 
state, avoiding crystallization, and to tune the degree of 
flexibility of the polymer chains, thus inhibiting complete 
folding of the chains. 

For the intermolecular interactions we choose 



~U{r,)j0{r,-r), (3) 

I 

ticular choice is done to obtain a polymer density in the 
melting state lower than the one achieved by consider- 
ing attractive interactions, thus avoiding crystallization 
at low temperatures j53| . e has energy units and controls 
the strength of the interaction, a has length units and 
fixes the typical length scale of the system. 
The intramolecular potential reads: 



if i, j n.n. 

U{rc) 0{rc ~ r), ifi,j'n.n.n. ^ ^ 

I 

ro/S. Miura et al. K5] found that the average cosine 
of the angle between two successive bond vectors along 
the chain equals 0.92, i.e., the chain is semi-rigid, for 
S ~ 2.5ro. 



B. The molecular motors 

A preliminary remark about the implementation of re- 
alistic molecular motors in a coarse-grained computer 
simulation study is obvious: no simple technique is able 
to take into account chemical activity in a molecular dy- 
namics calculation. Having this in mind, particular care 
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and physical insight has been put in the choice of a rea- 
sonable model for motor action on a particular monomer, 
say the ith in the ath chain, that we call ff^. Only 
a fixed fraction of polymers are provided with motors. 
This choice is mainly due to the requirements that mo- 
tor activity: i) should be homogeneously distributed in 
the sample; ii) should have a gentle but still detectable 
out-of-equilibrium effect on the behavior of the system; 
iii) should stay in the linear response regime. 

Following symmetry considerations, we chose to local- 
ize the non-conservative polymer activity at the core of 
the filament, i.e., we motorized the very central monomer 
in the polymer only. Note that, in general, this position 
does not coincide with the center of mass of the fila- 
ment and, most importantly, it is fixed. This situation is 
clearly at variance with the case of real systems, where 
molecular motors are able to slide along the filaments [6| . 
Nevertheless, we do not consider this a serious limitation 
to our approach, since we follow the dynamical behavior 
of the system on limited time scales over which motor 
activity can be indeed taken as localized. Also, we av- 
erage our results over a substantial number of different 
instances of motor realizations, mimicking in an indirect 
way motor motion relative to the filaments. A typical 
isolated flexible polymer configuration studied is shown 
in Fig. [H 




FIG. 1: The semi- flexible polymer model considered in this 
paper. It is formed by 21 monomers among which 20 are non- 
motorized (red spheres). The entire chemical motor activity 
is localized on the central monomer (blue sphere). Although 
all monomers have the same size we draw a larger central 
monomer in the sketch to emphasize its particular role. 

Similar to the case of self-propelled spherical parti- 
cles [i^, we mimic the motor effect with a time series of 
random isotropic kicks, generated by a suitable stochas- 



tic process controlled by three parameters: i) the num- 
ber of motorized polymers; ii) the intensity and direction 
of the motors forces; and iii) the activation time scale. 
During t steps of the molecular dynamics trajectory, in- 
dependent forces are applied to a fraction of randomly 
chosen motorized central monomers. The strength of the 
force exerted, = jfj^j is the same for all motors and 
is chosen to be a fraction of the mean conservative me- 
chanical force acting on the equivalent passive system, 
F = (NpN^r^ Eti Ifml- We quantify the mo- 

tor activity with the control parameter / — f^'^ /F. The 
direction of the motor force is chosen at random isotrop- 
ically, to avoid preferential flows in the system. The sub- 
set of propelled monomers and the directions of the ap- 
plied forces change at each power stroke of duration r. 
In short, we focus on adamant [s^ l motors, i.e. their 
action is completely independent of the structural rear- 
rangements they induce. 

C. Simulation details 

In a nutshell, given a model for the interaction forces 
between the different coarse-grained units forming a 
physical system, the technique of molecular dynamics 
simulation [s^ provide effective methods to numerically 
integrate the equation of motion, conforming to differ- 
ent statistical ensembles and macroscopic external con- 
ditions. It has been utilized for some decades now, to 
explore structural and dynamical properties of complex 
systems and it has been shown to perform quite success- 
fully also in the case of polymer systems [H, [H, 113] of 
interest here. 

In what follows we express all quantities in dimension- 
less reduced units: length, e nerg y and time-units are a, e 
and (mtT^/e)2, respectively [55|- AH the monomers have 
the same mass m — 1. In order to give an idea of typi- 
cal realistic orders of magnitude, we can associate these 
dimensionless quantities to reasonable values for the ex- 
perimentally measured counterparts [11] . We consider as 
a typical unit of energy, 2kBT. The typical lengths are 
measured in units of 0.4 nm, which gives a typical order 
of magnitude for forces of 20 pN at ambient temperature 
T = 300 K. 

Concretely, we integrated Eqs. (H]) numerically using 
Ermak's algorithm (ssl , [soj . with an integration time 
step 5t — 0.001 and a friction coefficient ^ = 10. This 
technique produces realistic configurations of the system, 
characterized by probability distributions conforming to 
the canonic NVT ensemble. We calculated statistical av- 
erages over several independent instances of the dynam- 
ics and we indicate them with angular brackets, (•.•). 
We used a cubic simulation box of total volume V with 
periodic boundary conditions. The parameters in the in- 
teraction potentials are: k = 9000 (the elastic constant). 
To = 0.4 (the equilibrium spring length) and S = 1. We 
simulated Np chains formed by Nm monomers each. The 
total number of interaction sites is = A^,, x Nm and 
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the number density p is given by the ratio N/V. The 
actual values for iV„ and Np have been fixed by some 
preliminary runs on the passive system detailed in the 
next section. 



III. PASSIVE SYSTEM: TUNING OF 
FILAMENT LENGTHS 

We performed some preliminary simulations to select 
the best suited values of the remaining parameters, Np 
and N,n in the limit in which the motors are switched off 
(passive system). We observe that the total number of 
interaction sites, N , is limited by our present computer 
capabilities, mainly due to the need of vast statistics. We 
first consider a system with Np = 250 polymers (a safe 
trade between the need for minimizing finite size effects 
and available computational resources) and values for iV,„ 
ranging from 5 to 101, to select the optimal polymer size. 

The wanted number density p=l is set by appropri- 
ately scaling the simulation box size. The target temper- 
ature, T — 0.8, corresponds to a safe liquid state. This 
choice of density and temperature has two important con- 
sequences. First, the passive system is in a liquid state 
characterized by relaxation times that are much shorter 
than the typical simulation total time scale and it is com- 
fortably in equilibrium with its environment. Therefore, 
averages of quantities that depend on only one time (after 
preparation) are time-independent and two-time correla- 
tions depend on the time-delay between the measurement 
of the two observables. More precisely, the dynamics are 
stationary. Second, the chosen thermodynamical point is 
well inside the stable liquid phase for the passive system. 
This allows us, in the active case, to safely assign any 
non-conventional structural or dynamical feature to the 
motor activity alone and not to some "normal" effect due 
to the interaction with some phase boundary. 



A. Collective static structure 

We start by quantifying the polymer size dependence 
of the structure of the system. In Fig. [5] we show the 
static structure factor 

5(q)=7V„liV-l(p(q,i)p*(q,i)), (5) 

with p{q, t) — X^ilTi Yl!l=i e*'' '''"'*-' the Fourier transform 
of the instantaneous density, r^a {t) is the time-dependent 
position of monomer i in polymer a. Isotropy implies 
S'(q) = S{q) and a spherical average over wave vectors 
of modulus q is considered. While the curve for Nm = 
5 clearly falls away from the remaining data, the static 
structure depends only weakly on Nm as shown by the 
collapse of all curves with Nm > 21. 

A detailed description of these results goes beyond the 
goal of this work, nevertheless we briefly discuss our find- 
ings in what follows. The S{q) of Fig. [2] show the typical 




Q 



FIG. 2: Study of the passive system at T = 0.8 and /3 = 1. 
The plot shows the collective static structure factor S{q) for 
systems with A'j, = 250 polymers with length Nm as indicated 
in the legend. For Nm larger than 21 no significant differences 
in the static structure on the entirely investigated wave- vector 
g-range are observed. 

features of a polymer melt, namely the presence of two 
principal peaks. The first one is usually related to inter- 
chains average distance and is at go — 4.8. The length 
associated to this value is about 1.26 which, as we will 
see below, is the typical distance between monomers per- 
taining to nearest neighbor polymer chains. The second 
peak, at position qi ~ 15.7 is associated to the equilib- 
rium bond distance tq. The feature of small intensity 
located at about g = 10 is not associated to any evi- 
dent length scale, being probably due to the convolution 
of density fluctuations of wave vectors between go and 
qi. Subsequent peaks, not shown in the figure, have less 
intensity reflecting the loss of spatial correlation in the 
liquid [60^. 

We note that the general shape of our curve is very 
similar to the one found in a united atoms model of 1,4- 
polybutadiene (see Ref. [5^ and references therein). In 
contrast, and as expected, S{q) is quite different from 
the one for highly flexible polymers modeled by a bond- 
spring model with Lennard- Jones interactions comple- 
mented by the finitely extensible non-linear elastic po- 
tential [6l|. 

/^From this preliminary analysis Nm = 21 appears as 
the minimal polymer length to be chosen to avoid any 
important dependence of the collective structure on the 
size of the polymer. 



B. Single polymer structure 

An equivalent characterization of the static structure 
of the passive system is based on the study of the pair dis- 
tribution function g{r) = V/n'^{YJ^^^ Hp^a ^('^ ~ '^^p)) 
with different choices of the sum ranges. This allows one 
to disentangle the inter and intra polymer correlations, as 



6 



shown in Fig.|3]for the case Np = 250 and Nm = 21. Data 
in panel a) are calculated summing over all monomers in 
the system; data in panels b) and c) , instead, correspond 
to the sums performed over monomers pertaining to dif- 
ferent polymers and monomers pertaining to the same 
chain, respectively. Data in a) are the convolution of the 
curves shown in b) and c). These results give additional 
support to the description proposed in Sect. IIII Al 
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FIG. 3: Study of the passive system at T = 0.8 and p = 1 
with Afp = 250 and Nra ~ 21. The pair distribution functions 
are computed summing over a) all beads in the system; b) 
beads pertaining to different polymers; c) beads pertaining to 
the same polymer. 

The analysis of the fluctuating radius of gyration [ioj 
averaged over all polymers 




confirms that the semi-flexible chains in the melt are non- 
Gaussian, in opposition to what is predicted by the stan- 
dard completely flexible polymer theory pol Isoj . This is 
shown in Fig. IHb), again for the case Nm = 21, where 
we plot the distribution of the radius of gyration for 
passive and active systems. The average value for the 
passive system, corresponding to Rg ~ 2.54 [Fig. [HI a)], 
has to be compared with the Gaussian model value, 
i?,, =r-oyAW6-0.75. 

C. Equilibrium dynamics 

The question of the role played by rigidity in the dy- 
namics of a concentrated melt of linear polymers is highly 
non-trivial (see, among others, Ref. [62]) and a compre- 
hensive study of this issue goes far beyond our present 
interests. Here we are primarily interested in a qualita- 
tive characterization of the self dynamics of the passive 
polymer system. In particular, we need to check if the 
value Nm = 21 which appears to be the best choice for 



a set of Np = 250 interacting filaments is characterized 
by a structural dynamics taking place on reasonable time 
scales. Our main results are summarized in Figs. |4] and [5l 
In Fig. 0] a) we show the averaged (three-dimensional) 
mean-squared displacements 

^'it) = -{J2T.\'^^a{t)-rU0)\'), (7) 
a—1 i—1 

and in panel b) we show the self-intermediate scattering 
functions 

-, Np N,„ 
a— 1 i—1 

at q = 4.8, the wave- vector of the first diffraction peak 
in the static structure factor, in both cases calculated 
considering all beads and as a function of t. The chain 
lengths are indicated in the figure. Our data are con- 
sistent, as expected, with a picture where the larger the 
polymer size, the slower the relaxation. We also note 
that these data show the typical signatures of a liquid 
state, with no dynamical arrest nor crystallization. 
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FIG. 4: a) Mean-square displacement and b) self intermediate 
scattering function in the passive system with p = 1 at T = 
0.8 as a function of t for different values of Nm given in the 
legend of panel b). In all cases Np — 250 and all monomers 
have been taken into account. The limited time range used 
still allows for the determination of a self diffusion coefficient 
from a linear fit of the long time part of the mean squared 
displacement. Some examples are shown with solid lines. The 
self intermediate correlation function is calculated at a value 
of g = 4.8 corresponding to the first diffraction peak of the 
static structure factor. 

The mean-square displacement in Fig. 2] a), presents 
the typical bending separating short time-scales {t ^ 

10^^) with ballistic motion from long time-scales (t ^ 
10^^) with, at least for short chains, diffusive behavior. 
Indeed, although the investigated time range is quite lim- 
ited, for N„i < N^ = 50 we can estimate a diffusion co- 
efficient from the Einstein relation D = A^(i)/6i in the 
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limit t OQ. These data are shown in Fig. [5] a), and 
points for < are consistent with a power law 
iV~^ (dashed line). Data for Nm > N*^ are much less 
convincing, saturating at a constant value. This is pos- 
sibly due to the limited time scale reached, which does 
not actually allow us to discriminate between a diffusive 
behavior and a Rouse-like behavior oc t^-^^ [1^, [s^l 
which should be a signature of limited entanglement. 
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FIG. 5: a) Self diffusion coefficient D calculated by a lin- 
ear fit of the long time linear behavior of the mean squared 
displacement, as explained in the text. The dashed line is 
a guide-to-the-eye. b) The relaxation time extracted from 
Fs{q,t) for g=4.8 as described in the text. Note the apparent 
dynamical crossover at ~ 50. We remind that all quanti- 
ties are calculated from all monomer coordinates, as discussed 
in the text. 

The quality of our data for the intermediate scattering 
functions [Fig. 0] b)] allows us to extract quite comfort- 
ably an estimate for the structural relaxation time by 
using the model independent definition Fs{q,Ta) = I/e. 
Our results are shown in Fig. [5] b) . We find evidence 
for a dynamical crossover similar to the Rouse-like to 
reptation-like one [50]. More precisely, we find a relax- 
ation time Ta — Nm^ for Nm < N*^ and ~ Nm^ 
for Njn > N*^, with ~ 50. We tentatively asso- 
ciate the anomalous power-law dependencies to the fact 
that we are dealing with semi-flexible polymers. The 
wave- vector dependence of Ta is the usual hydrodynamic 
behavior Ta — (not shown). 

We conclude our brief study of the chain-length de- 
pendence of the structure and dynamics of a passive 
semi-flexible polymer melt by stating that, surprisingly 
enough, to the best of our knowledge, this problem has 
not been fully investigated in the past. The iV,~^ depen- 
dence of the self-diffusion coefficient is a known scaling 
law of quite generic polymer dynamics but we know of no 
previous numerical studies to compare with concerning 
the Nm dependence of r^. A complete characterization 
of a possible dynamical crossover corresponding to a par- 
ticular chain size is also lacking. All these matters are of 
limited interest in the context of the study of Tgy f but 
call for further work. Here we conclude by fixing the final 



parameter Nm = 21 corresponding to a passive system 
with quite standard structural and dynamical features, 
taking place on time scales comfortably reachable in nu- 
merical experiments. 



IV. ACTIVE SYSTEM: THE EFFECT OF 
MOTOR ACTIVITY 

We now turn to the investigation of the structure 
and dynamics of the active semi-fiexible polymer melt 
in which motors are switched on. We stress here that 
external macroscopic conditions are fixed and the only 
control parameter is the motor activity. We recall that 
we quantify the motor activity via a single parameter 
/ = f^^ /F, i.e., the ratio (in the range 0.1 -;- 1) between 
the force exerted by the motor and the mean mechan- 
ical force, F ~ 163.5, acting on the equivalent passive 
system. We fix the activation motor time as r = 500 
molecular dynamics steps. Similar to the inactive case 
detailed in Section IIII[ we consider Np — 250 polymers 
with Njn = 21 monomers each, for a total of A^ = 5250 
interaction sites. 




FIG. 6: Typical snapshots of the system of motorized fila- 
ments at the indicated values of the motor activity, /. Brown 
beads indicate the instantaneous positions of the molecular 
motors. 

In Fig.[6]we show typical snapshots of the system pro- 
duced following the above prescriptions, at the indicated 
values of motor activity. Differences are not immedi- 
ately evident and we need to characterize quantitatively 
the effect of motors on intramolecular and intermolecular 
structure. 
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A. Thermodynamics 

We start by analyzing the /-dependence of a few ther- 
modynamic observables. Figure [7] presents the / de- 
pendence of the inter-molecular and intra-molecular en- 
ergy, see Eqs. ^ and respectively, as well as the 
pressure defined through the usual virial relation PV — 
NT+l/3{J2ia^ia'^ia) ^M- The thrcc quantities incrcasc 
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FIG. 7: Thermodynamical observables for the active system, 
as a function of the motor activity / defined in the text. From 
left to right we show the inter-molecular repulsive Lennard- 
Jones contribution a), the (total) intra-molecular energy b) 
and the total pressure c) (defined in Eqs. (Ol and Q and 
calculated from the usual virial relation, respectively). 

smoothly for moderate forces, / ^ 0.2, next they cross 
over to an approximately linear, in the case of Umter, 
and sigmoidal, in the cases of Umtra and P, shape. We 
note that the overall variation of the intermolecular po- 
tential energy is of order 8, while it is of order 0.1 for 
intramolecular (sum of the repulsive Lennard Jones and 
bonding) energy and total pressure. We also note that for 
the latter only the configurational term plays a role, be- 
ing T constant. The results of panel a) seem to suggest 
a picture where distances between beads pertaining to 
different filaments decrease with increasing /, therefore 
indicating crowding. 

The pressure results are also very interesting. Indeed, 
pressure is related to the (negative) diagonal components 
of the stress tensor. Thus we can anticipate that motor 
activity has a non-trivial effect on the mechanical prop- 
erties of the filaments. This is actually what is observed 
in real systems, see e.g. Ref. [63[. Already at the level 
of simple thermodynamic quantities related to configura- 
tional properties motor activity is non negligible and we 
expect to find sensible implications on the static struc- 
ture as well. 



B. Collective structure 

In this section we consider the static structure factor 
S{q), see Eq. ([S|), and the gyration ratio defined in Eq. ^ 



and we try to give a consistent picture of the motor ac- 
tivity effect on the static structure of the system. 

In Fig. [5] we show the static structure factor as a func- 
tion of motor activity. All curves have been shifted by 
a fixed amount for clarity, / increases from bottom to 
top. The general features of the structure factor are not 
altered by the action of the motors but, upon careful 
inspection, one observes some subtle changes. 

Upon increasing the strength of the motor activity, the 
main peak position shifts to higher values of q, meaning 
that the average nearest neighbor distance between beads 
pertaining to different polymers decreases (recall that the 
first and highest peak is mainly determined by interchain 
correlations). This is also consistent with the crowding 
picture suggested from the study of the thermodynamic 
quantities in Sect. IIV Al The width of the first peak in- 
creases while its height decreases with increasing /. This 
trend is similar to the one that one finds when increasing 
the bath temperature in a passive polymer liquid, see for 
instance the simulations in [6lj. The motors have then 
the mixed effect of making the melt more compact and 
more disordered simultaneously. 
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FIG. 8: Static structure factor in the active sample at the 
investigated motor activities parametrized by / (highest / on 
top). Data have been shifted vertically for clarity. The first 
diffraction peak increases with / as described in the text. 

The second and third peaks in S{q) follow the same 
trend, see Fig. HI suggesting that the effect of motor ac- 
tivity is to fold the chains. This is consistent with the 
fact that the radius of gyration decreases with increasing 
/, see Fig. Ha). 

Working in real space one corroborates this picture. 
Not much information can be extracted by simply looking 
at the system snapshots of Fig. [51 In particular, the fact 
that increasing / has a sensible influence on the gyration 
ratio of the filaments is not evident from the snapshots, 
as it is when we quantify the analysis in Fig. [HI We show 
in Fig. [3] a) the average value of the gyration ratio as a 
function of motor activity. It decreases from 2.54 for the 
passive system to 1.7 for / = 1. 

The probability distribution of i?g [Fig.[S]b)] suggests 
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FIG. 9: a) Averaged gyration ratio, b) Probability distribu- 
tions of the gyration radius at all investigated motor activi- 
ties. The structure is highly non-trivial with an average value 
shifting towards lower values for increasing /. A more precise 
characterization of the distribution is given by its variance c), 
skewness d) and kurtosis e). 



an even more complex scenario. In the passive system 
it is non-Gaussian and very asymmetric. We can char- 
acterize these curves v ery precisely with the variance, 
skewness and kurtosis [gj] dependence on /, as shown 
in Fig. [H]c)-e). A cross-over between regions of different 
curvature is evident, but we are not in the position to 
justify this complex behavior. It would be very inter- 
esting to compare these data to analogous experimental 
measurements using scattering techniques. 

In conclusion, motor activity apparently has the effect 
of pushing closer the filaments, which at the same time 
fold substantially. There is a competition between two 
effects: crowding, which is related to excluded volume ef- 
fects and tends to slow down the dynamics; and folding of 
the chains, which dynamically lifts the topological con- 
straints, therefore decreasing entanglement. This must 
have important consequences on the long-time evolution 
of our system, as we will see in the next section. 



C. Dynamics 

In this section we quantify the effect of motor activity 
on the dynamics and we compare it to the relaxation 
in the passive limit. We stress again that we work at 
constant volume and temperature, which means that all 
observed differences in the dynamical behavior are only 
due to the motor effects which are in any case chosen to 
be mild. 

Figures [10] and [TT] (main panel) display the mean- 
square displacement, Eq. ([7]), and the self-intermediate 
scattering function, Eq. (|5]) - for the wave- vector qq ~ 4.8 
of the main peak in S{q) in the passive case - as a func- 
tion of time-delay for all investigated motor forces given 
in the key. The time range considered is quite extended 
and allows us to comfortably reach the diffusive regime 
of the global displacement in which the self-intermediate 
scattering functions safely decay to zero. 
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FIG. 10: Mean-squared displacement in the active system 
with the motor activities given in the key (higher activities 
on top). 

Qualitatively, the data resemble the ones found in the 
passive system, see Fig. |4l but depend now on the mo- 
tor strength. All the data are compatible with a picture 
where the collective dynamics of filaments gets faster un- 
der stronger motors. This is also consistent with our 
results on the gyration ratio and supports our interpre- 
tation: the augmented folding of the filaments seems to 
solve local topological constraints and diminishes the ef- 
fect of entanglement. 

Additional interesting information is included in the 
data of Fig. [TTJ Solid lines are best fit of our data to a 
curve of the form exp— (i/r)^. Values for the (3 coeffi- 
cient as a function of / are included in the inset. Data in- 
creases with /, implying that the relaxation crosses over 
from stretched exponential (/? < 1) to compressed expo- 
nential {(3 > 1). This unusual shape for the intermedi- 
ate scattering function has been already observed in soft 
matter systems, like colloidal gels, emulsions or Laponite 
suspensions, as surveyed in a number of reviews, see, e.g., 
Ref. [1^. In these references, the cross over is related to 
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FIG. 11: Main panel: Time-dependence of the intermediate 
scattering functions at the wave-vector corresponding to the 
first maximum of the static structure factor at the investi- 
gated motor activities (activities increase to the left). Inset: 
values of the exponential parameter beta estimated from a fit 
of the data to the form exp— (t/r)'' (see text for details). 
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FIG. 12: Diffusion coefficient D and structural (self- 
)relaxation time r as a function of /. The relaxation time 
values have been estimated from a fit of the data to the form 
exp — (t/r)*^ (see text for details). The data follow power laws 
with exponents indicated in the figure. 



randomly distributed internal stress sources acting on the 
sample. It is tempting to associate this sources to motor 
activity but, again, it is difficult to address this subject 
here. 

An even more quantitative analysis is presented in 
Fig. [T2 in which the dependence of the diffusion coef- 
ficient and the (inverse) relaxation time (estimated by 
the fit detailed above) on the strength of the motors is 
analyzed. Data are shown renormalized to the value of 
the passive system. It is very interesting to note that 
both sets of data are consistent with two very simple 
functional forms (solid lines in the figure), D/Df=o ~ 
1 + 1423 X /2-29 and (r/r/^o)"^ ~ 1 -t- 19 x /, respec- 
tively. Note the extremely large variation of the diffusion 
coefficient that extends over three decades, and the less 
pronounced variation of the inverse relaxation time. 

An observation is in order at this point. The /- 
dependence of the diffusion coefficient is reminiscent of 
the findings discussed in Refs. [4l| and [i^. For the ac- 
tive colloidal particles used in these studies the effective 
diffusion coefficient can be written as D ^ Do + l/3 v^tji, 
where Do = Df=o = ksT / {Q-K-qR) is the diffusion coeffi- 
cient of the passive system, v is the velocity of the par- 
ticle and Tji = SirrjR^ /{ksT) is the rotational relaxation 
time. R and 77 are the typical radius of the colloidal par- 
ticle and the viscosity of the solvent, respectively. This 
relation can also be written as D/Dq = (1 + l/S^'e^), 
where Pe is the Peclet number, defined as Pe = vR/Do 
and characterizes the particle activity [4l|, |63| • 

In our case the exponent (2.3 ± 0.1) is close to 2. The 
difference could be ascribed to the uncertainties in the 
determination of D or to the presence of corrections to 
the correct scaling. It is tempting to imagine that / is in 
this case playing the role of the Peclet number. We will 



see in our final discussion that this idea is indeed rea- 
sonable and our results on the effective temperature and 
analogous experimental data j4l| seem to be compatible 
with the conjecture. 

We have not been able to find a rationale for the linear 
dependence of the inverse of the relaxation time, which 
is also a quite striking result. 



V. EFFECTIVE TEMPERATURE(S) IN 
COMPLEX ACTIVE MATTER 

The analysis above proved that both structure and dy- 
namics of the active system are influenced by motor ac- 
tivity in a very complex fashion. We ask now whether it 
is possible to embed all this complexity in a single param- 
eter, also prone to be directly determined in experiments. 
The chosen parameter is the effective temperature T^j / a 
quantity already well characterized in different contexts, 
as we have discussed in the introduction [30| . 

In equilibrium, a model independent relation exists be- 
tween the linear response of an observable, say A, to 
a perturbation applied to the system and the correla- 
tion between the same observable and the one that "re- 
ceives" the perturbation, say B. This is the fluctuation- 
dissipation theorem (FDT) that reads 

XAB{ti,t2) = ^[CAB{tl,h)~CAB{tl,t2)], (9) 



where 



XAB{tl,t2) 



ShBit') 



hB=0 



CABihM) - {A{h)B{t2)), 



(10) 
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FIG. 13: Fluctuation-Dissipation relations in the motorized 
polymer melt. The red data on top correspond to the passive 
system. Solid lines are linear fits used to estimate the value 
of the effective temperatures. The dashed (red) line has slope 
— l/T, thus confirming the validity of the FDT at equilibrium 
for the passive system. 



and the perturbation is such that H ^ H — hsB. 
The prefactor is the inverse of the temperature 

of the bath with which the system is in contact. In 
equilibrium CyiB(ii,ti) is a constant, both XAB{tiTt2) 
and CAB{ti,t2) are just functions of the time difference 
t ^ ti — t2, and Eq. ^ becomes 

XABit) = ^[C'ABiO) - CABit)] , (11) 

for i > 0. In the frequency domain the FDT reads 
2T'^Rab{i^) — vCab{i^) with the sign convention for the 
Fourier transform such that RABiy) — dt e^"*^ RAB{t)- 



another one characterizing the structural relaxation at 
long time-delays. 

In interacting particle problems [69| it is custom- 
ary to use A{t) = l/iVX;.^iQe*''-^'(*), and B{t) = 
2 X^iLi cos[(f • fiit)], where the field = ±1 with prob- 
ability a half. With this choice the relevant equilibrium 
correlation function is the intermediate scattering func- 
tion of Eq. ^ once the results are averaged over many 
instances of the e^. We have also used these observables 
with q — qo and we averaged over 10^ independent field 
realizations. 

We performed these measurements for different values 
of the forcing / and the results are shown in Fig. [T^] in 
the form of the parametric plot x{t) against C{t). The 
curves clearly show two regimes with a relatively smooth 
crossover between them. We do not attribute particular 
importance to the first regime, corresponding to small t 
or high frequencies, which represents the fast vibration of 
monomers, surrounded by their environment, before any 
structural relation takes place. The second regime, being 
the one of long t and small frequencies, is more interesting 
since it describes the actual structural relaxation in the 
sample. 

In all cases, the parametric plot in this regime is rather 
well described by a straight line [72J, implying that T^ff 
is a constant over this time-scale. In the figure the nu- 
merical data are represented with points and the lines 
are the result of fits to linear functions performed on the 
long-times part of the curves. For the passive system 
the construction yields a straight line with slope equal 
to —l/T, confirming that the FDT holds in this case. 
The active systems have decreasing absolute value of the 
slope for increasing forcing and will be compared to other 
measurements of Te// below. Data from the x(C) deter- 
mination of Teff are reported in Fig. 1201 with open (red) 
squares where we also confront to other measurements as 
we will discuss in depth later. 



A. Correlation-response parametric plot 

One way of measuring the effective temperature in a 
non-equilibrium steady state consists in using the para- 
metric relation between the out-of-equilibrium integrated 
linear response, XAB{t), and the correlation of a chosen 
pair of observables, C^B(t), and associating minus its 
inverse slope with a possibly time-dependent parameter 

Teffit) mM- 

XABit) = TfAjJ-JCABiO) - CAB(t)\ ■ (12) 

Such measurements were performed in different aging, 
that is to say non-stationary, glassy systems (for details 
see ^) and in out of equilibrium steady states (see, 
among others, Refs. |69l - [7l| ) demonstrating that Te//(t) 
is piecewise constant taking, typically, two values: the 
one of the environmental bath at short time-delays and 



B. Tracer particles 

The construction detailed in the previous section is 
very powerful and provides us with direct unambiguous 
results. Indeed, it is at the same time a tool to demon- 
strate the existence of an effective temperature and an 
operative method to calculate it (69j . The method is eas- 
ily exploitable in computer simulations, where the micro- 
scopic implementation of the construction does not pose 
particular problems. The situation is completely differ- 
ent in experiments, where the coupling of each (complex) 
unit forming the system to the perturbation field is far 
from obvious. This difficulty has been circumvented by 
considering the dynamics of tracers (or intruders) of sim- 
ple geometry, e.g., spheres. These objects of micrometric 
size are nowadays quite simple to control and manipulate 
by active microrheology techniques. When immersed in 
extended complex environments they have been demon- 
strated to couple to them as to an external reservoir. 
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Their dynamics provide a wealth of information related 
to the underlying matrix state. Below we discuss a few 
methods to implement these ideas in our simulations, 
based on the use of free and driven tracer particles, re- 
spectively. A typical snapshot of our filament system 
with embedded tracer particles is shown in Fig. [141 We 
have used spherical particles that, for visualization conve- 
nience, we draw with an augmented radius in the figure. 




FIG. 14: Molecular dynamics computer simulation imple- 
mentation of the experiment with tracer particles. We show 
a typical polymer configuration with 10 tracer particles; 50 
per cent of polymers are motorized. Tracers have been drawn 
with an augmented radius for clarity. 



1. Free tracers 

The effective temperature should be measurable with 
a fine tuned thermometer that couples to the structural 
relaxation of the sample. A possible realization of an 
adapted thermometer is constituted by a massive - i.e., 
of mass much larger than the mass of the polymer beads ~ 
free tracer particle immersed in the system. Indeed, it is 
possible to show that the kinetic energy of the tracer can 
be related to the system's effective temperature via an 
equipartition theorem for slow modes. This kind of mea- 
surement were suggested in p9j and they were performed 
experimentally in granular matter [ssj . colloidal suspen- 
sio ns i40.] and numerically in models of granular mat- 
ter [71| and sheared atomic Lennard- Jones mixtures [1^ . 

We consider Ntr = 10 tracer spherical particles of vari- 
able mass rritr in the interval [10 — 10^], which inter- 
act with all monomers in the sample via the repulsive 
Lennard- Jones potential of Eq. Tracers are not cou- 



pled to any external thermal bath and their dynamical 
evolution is controlled by the usual Newton equations of 
motion that we have integrated with the velocity- Verlet 
method with St = 0.001 in the NVE ensemble. The dy- 
namics of polymer beads is still treated by a Langevin 
approach. Note that the tracers do not interact among 
themselves and the use of Ntr tracers allows for better 
statistics. 

The results of our measurements for / = 0.5 are shown 
in Fig. [ini Panel a) displays the probability distributions 
of the tracer velocity, where d = |v| and its depen- 
dence on the tracer's mass. The data are shown using the 
scaled variables p{v) / mll"^ against vrriH'^ . The distribu- 
tions are very close to Gaussian for all masses. The data 
are shown with points and the lines are fits to this func- 
tional form from which we extract an estimate for the 
effective temperature as a function of the tracer's mass 
that we show in panel b). 
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FIG. 15: a) Rescaled probability distribution of the veloc- 
ities of the massive tracer particles immersed in an active 
sample with / — 0.5. The distributions are remarkably Gaus- 
sian for all tracer masses. Their variance increases with in- 
creasing mass, b) The effective temperature, extracted from 
a Maxwell- like fit to the distributions shown in panel a) , as a 
function of the tracers' mass mt,-. 

The interpretation of these results is quite clear. Light 
tracers quickly react to the rapid motion of the polymers 
and test the short time-delay dynamics of the sample. In 
contrast, heavy tracers are insensitive to such rapid mo- 
tion and displace themselves only due to the large-scale 
structural reorganization of the sample, that is to say, the 
long time-delay dynamics of the embedding matrix. The 
interesting long times - or low frequency - value of the 
effective temperature is then accessed in the large mass 
limit of this measurement. 

^From the figure we see that at mtr — fO'^ the varia- 
tion of Teff stops and the curve saturates to a constant 
level from which we read Teff{f = 0.5) ~ 5.1. Simi- 
lar measurements for other values of the motors strength 
yield the data collected in Fig. [201 with open (blue) cir- 
cles. In what follows we will consider mtr = lO'^ which 
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FIG. 16: a) Uni-dimensional mean-squared displacement of 
the free tracer particles at different values of the motor activ- 
ity, b) Parametric plot of the tracer particles mean squared 
displacement as a function of the driven displacement, at the 
indicated values of motor activity. The effective temperatures 
have been estimated with a linear fit (solid lines) to the large 
mobility part of the curves. 



corresponds to the smallest tracers mass coupling to the 
relevant slow modes of the filaments. 



2. Driven tracers 

Yet another independent measurement of T^f / is based 
on the use of driven tracer particles. More precisely, it is 
given by the study of diffusion and mobility of an ensem- 
ble of non-interacting tracer particles immersed in the 
active sample. These two quantities are related by the 
Einstein relation, one particular form of the fluctuation 
dissipation theorem. The (free) mean squared displace- 
ment, is the relevant equilibrium correlation function; the 
associated response function is the displacement induced 
by applying a small constant external force of intensity 
h to the tracer. Due to isotropy, we focus on the one di- 
mensional version of these observables and consider the 
external force oriented in the x direction. 

In Fig. [in] a) we show the one dimensional free 
mean-square displacement for tracer particles, A^(t) = 



l/^tr(E!^I(2^r(0 - xf{Q)f), in the x direction as a 
function of time-lag, for different forcing strengths given 
in the key. Results in the other two directions are equiv- 
alent. In all cases the tracer particles mean-square dis- 
placements cross over from ballistic motion at short time- 
delays to diffusive motion at longer time-delays. These 
two limits are shown with straight lines as guides-to-the- 
eye in the figure. 

The one-dimensional mobility in the x direction is 
calculated as Xx{i) = ^x{t)/{ht), where Ax{t) = 

i /Ntr{J2f=ii^T i^) ~ ^Ti^))) is the displacement calcu- 
lated in independent simulations by applying to the trac- 
ers a small force h = O.lFx ~ 16. Si. 

Conforming to the Einstein relation D^/xx = T, with 
Dc the diffusion coefficient in the x direction and Xx the 
saturated mobility, the effective temperature should also 
be the parameter linking the tracers' mean-square dis- 
placement and their driven displacement under a weak 
external applied force: 



eff- 



Axit) 



(13) 



in the long time-lag limit in which the dynamics is dif- 
fusive. This equation is plotted in parametric form 
in Fig. [16] b) for the passive system and seven motor 
strengths given in the key. The effective temperature is 
given by the slope of these curves, as shown in the fig- 
ure. These data are included in Fig. [5D] as open (green) 
triangles. 
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FIG. 17: The effective temperature, defined as the ratio 
/iA^(!/)/(2Aa;(!/)), as a function of frequency. We represent 
as error bars the imaginary part of the ratio, to quantify the 
precision of the numerical Fourier transformation. Diamonds 
on the left are the data extracted from the linear fit shown in 
Fig. [16] b). 
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3. Frequency- dependent Tt 



Parametric plots in the time-domain (Figs. [T3landfT6l 
b)) are the standard way of identifying the values taken 
by the effective temperature in well-separated dynamic 
regimes each with a constant value of T^ff. They may 
hide, though, long crossovers in which T^^ff does indeed 
depend on time. Besides, the determination of an ef- 
fective temperature in experiments very often relies on 
measurements of spontaneous and induced fluctuations 
in the frequency domain. 

In Fig. [17] we represent the same data shown in 
Fig. Iinib) in a closer way to what is implemented in ex- 
periments. We have Fourier transformed both the trac- 
ers' mean-square displacement and displacement under 
an external field and we present the ratio. 



(14) 



as a function of frequency. This determination of T^ff 
should coincide with the one of the time-domain in the 
regimes in which Te// is constant. And this is indeed 
what we have found: for all the activation strengths and 
in the accessible frequency range, the curves show a rel- 
atively constant force-dependent plateau. In the limit of 
very low frequency we recover the frequency-independent 
values (diamonds) obtained by the linear fit of the para- 
metric plot in the time-domain (Fig. [T6|) . Although the 
accessible frequency range is limited by the sampling rate 
(upper bound) and the total time length (lower bound) 
of the trajectory and in principle one cannot exclude a 
different behavior outside this region, we conclude that 
no evidence of frequency dependent effective temperature 
has been found in the present system. 



4- Free tracer dynamics. 
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FIG. 18: Parametric plot of the rescaled diffusion coefficients 
for a free tracer at different values of the number density p 
and bath temperature T as a function of the rescaled effective 
temperatures. The dashed line is a linear fit to the data. More 
details are given in the text. 



the effective temperature calculated by the massive par- 
ticle method. This means that once the master curve is 
known, we can predict the value of T^f / by just following 
the free dynamics of the tracers {i.e., by computing D). 

We conclude this section with the following observa- 
tion [tI]. In the case of out-of-equilibrium glassy systems, 
one-time quantities converge to constants and, in par- 
ticular, thermodynamic observables receive contributions 
from the fast regime at the ambient temperature and the 
slow regime at the (higher) value taken by the effective 
temperature, both as if they were equilibrium ones. In a 
sense, for these observables, the out-of-equilibrium sys- 
tem behaves as if it were equilibrated at T^ff- We could 



The power of methods based on intruders has been 
largely demonstrated in the literature [t^ and in the 
present work. The results contained in Fig. [TH] confirm 
this fact and also allow one to add an interesting infor- 
mation: knowledge of the effective temperature depen- 
dence on / in particular macroscopic external conditions, 
i.e., temperature or pressure (density), allows one to pre- 
dict Tf.f f at any different macroscopic external conditions 
with one free tracer measure only. Indeed, we have found 
that a simple power law relates the diffusion coefficient 
of a free tracer, D, to the effective temperature charac- 
terizing the polymer matrix. All data corresponding to 
different temperatures or densities can be collapsed on a 
master curve by a simple rescaling: 



Dp 
T 



T, 



T 



(15) 



with Of ~ 1. This is exhibited in Fig. [T5] where we 
plot the diffusion coefficients of tracer particles against 
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FIG. 19: Comparison between gyration radius of passive and 
active systems at the indicated temperatures and motor ac- 
tivities. See discussion in the text for details. 
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therefore ask, what can we say about the structure of 
the passive system when equihbrated at T^ffl In Fig.fTOl 
we plot the probabihty distributions of the radius of gy- 
ration (a one-time quantity), at the indicated values of 
bath temperature and motor activity. It is evident that 
the distributions at the higher bath temperatures for the 
passive case (open symbols) nicely reproduce those at 
T = 0.8 and the indicated values of activity (closed sym- 
bols). The problem is that these bath temperatures are 
not those we could predict from the data in Fig. [501 We 
believe that the connection between effective tempera- 
ture and static structure is even more problematic that 
one could expect. 



tend to Teff = T in the limit of vanishing activation, 
as they should, and the deviation from the environmen- 
tal temperature monotonically increases with increasing 
forcing. The data can be fitted by the empirical law 
Teff/T = 1 + with 7 = 15.41 for filaments and 
7 = 1.18 for particles, showed with dashed lines in the fig- 
ure. We note that this finding seems to support our con- 
jecture that the parameter / plays here a role analogous 
to the Peclet number for colloidal active particles used in 
a few experiments [4l|, as discussed in Section lIV CI 



VI. CONCLUSIONS AND PERSPECTIVES 



C. Survey on Te// vs motor activity. 

In conclusion, in this section we have shown that T^f / 
characterizes the chemical motor activity which we have 
schematized by the single parameter /. Figure [20] col- 
lects all our results both for self-propelled particles [43] 
and motorized semi-flexible filaments. The data are pre- 
sented in the form T^f f/T against / and we also show (di- 
amonds) average values with error bars calculated from 
the different determinations of T^f f . 
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FIG. 20: This figure summarizes the entire set of data that 
we have generated, both for spherical self-propelled parti- 
cles [i^l and motorized semi-flexible polymers, on a linear- 
logarithmic scale. The different sets of data are indicated in 
the legend. Empirically, we find that all data are compatible 
with a power-law dependence of the effective temperature, 
Teff/T ^1 + -yf (dashed lines). 

The figure demonstrates that independent determina- 
tions of Teff yield consistent results suggesting that, in 
these systems, there is a T^ff parameter with a ther- 
modynamic meaning. The data for semi-flexible lines 
and spherical particles fall on different curves but there 
was no reason for universality in this sense. They both 



The limits of the approach used in the present work are 
clear. The model we have developed is quite crude and 
does not take into account explicitly any chemical activ- 
ity nor particularly complex molecular structure. More- 
over, we have chosen the most favorable conditions for 
producing and observing the desired effects, at the risk of 
perhaps being slightly unrealistic. Nonetheless, we think 
that our results are relevant for understanding very com- 
plex phenomena in real systems. This conviction is based 
on two facts. 

First, our work constitutes a genuine microscopic ap- 
proach to active matter. In particular, no uncontrolled 
hypothesis has been made, the physics of the passive sys- 
tem being under complete control: it conforms to previ- 
ously known results and it is only due to the effect of con- 
servative forces of very general nature. Also, the imple- 
mentation of the molecular motors, the exclusive source 
of non equilibrium fluctuations in the system, is suffi- 
ciently general to be representative of, at least, adamant 
motors. In conclusion, we expect our results to be quite 
universal. Second, as we discuss below, our results are 
reminiscent of experimental findings reported in the lit- 
erature. From this point of view the approach presented 
here can help both in rationalizing experimental obser- 
vations and in suggesting new research directions. 



A. Context 

A number of recent publications explored fluctuation 
dissipation relations in biological systems or biologically 
motivated models. We devote the following paragraphs 
to the discussion of results contained in the papers that 
are relevant in the context of the present work. 

Ziebert and Aranson Wd\ performed Langevin simula- 
tions of two-dimensional solutions of semi-flexible polar 
filaments interacting via molecular motors. They stud- 
ied the rheology and structural properties in diluted sam- 
ples and identified an "active temperature" Ta via a phe- 
nomenological scaling of the elastic and loss moduli. Ta 
increases with the average number of motors per filament 
and with the motor velocity. This definition is, though, 
not necessarily equivalent to the one that we used. Closer 
the our work are the following studies. 
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By comparing a hair bundle's spontaneous oscillations 
with its response to small mechanical stimuli, Martin 
et al. [tHI demonstrated the breakdown of the FDT, 
and thus confirmed that a hair bundle's spontaneous 
movements are produced by energy-consuming elements 
within the hair cell. 

In [63j the fiuctuation-dissipation relation in a model 
cell system formed by three fundamental components 
(myosin II, actin filaments, and cross-linkers) was studied 
in equilibrium and out-of-equilibrium conditions. FDT 
holds in the former and fails in the latter situation. In 
these papers, the violation of FDT was used as a proof of 
the non-equilibrium character of the dynamics induced 
by molecular motors. 

The effective temperature stemming from deviations 
from FDT has been used to study the stability of mo- 
torized particles in cytoplasm jsj] ■ Interestingly enough, 
values of Tg/ / that are lower than the ambient tempera- 
ture were found when considering motors that are aware 
of the free-energy landscape and act only in the direction 
of lowering its instantaneous value. The motors used in 
our study are not of this kind and the values of T^f / ob- 
tained are consistently higher than the ambient T. Such 
an "inversion" also occurs in cases in which the initial 
configuration is chosen to be to be one of equilibrium at 
a lower T than the working one j77| . 

A simple birth-death process and a general two-species 
interacting process, as simplified models of gene net- 
works, were studied from the effective temperature per- 
spective in [78| . The aim of this paper was to propose to 
use Tef / as a means to distinguish between intrinsic and 
extrinsic noise sources. 

Morozov et al. [t^ studied a model of the cytoskele- 
tal network made of semi-flexible polymers subject to 
thermal and motor-induced fluctuations, modeled as ei- 
ther additive or multiplicative colored noise. They found 
violations of the conventional FDT leading to a T^ff 
that exceeds the thermodynamic temperature T only in 
the low-frequency domain where motor agitation prevails 
over thermal fluctuations. This is consistent with our 
numerical findings. However, T^ff in this phenomeno- 
logical model is continuously frequency dependent while 
our data suggest that there is a low frequency regime in 
which it reaches a constant value, as in glassy systems. 

Fluctuation-dissipation ratios were used to extract in- 
formation about the degree of frustration, due to the ex- 
istence of many metastable disordered states, in two self- 
assembly processes: the formation of viral capsids, and 
the crystallization of sticky discs [SO^ . 

Joly et al. [42] used numerical techniques to study 
the non-equilibrium steady state dynamics of a heated 
crystalline nanoparticle suspended in a fluid. A two- 
temperature scenario holds in this case. Indeed, by com- 
paring the mobility to the velocity correlation function, 
FDT approximately holds at short-time lags with a tem- 
perature value that coincides with the kinetic one. In 
contrast, at long-time lags data are compatible with the 
temperature estimated by using the Einstein relation. 



Ben-Isaac et al. [8l| analyzed the Langevin dynamics 
of a single particle randomly kicked by motors as a toy 
model of active matter. They focused on the deviations 
of Teff from T and the distribution of velocities from 
Gaussian as parameters that quantify the out of equi- 
librium nature of the dynamics, flnding that while the 
former is robust the latter is not. Comparison to mea- 
surements in ATP-driven red-blood cells is also provided 
in this paper. 

Finally, let us sum up the results in the experimental 
work of Palacci et al. [41| that is, indeed, the closest to 
our work. These authors investigated the effective tem- 
perature by following Perrin's analysis of the density pro- 
file in the steady state of an active colloidal suspension 
under gravity. The active particles used - Janus parti- 
cles - are chemically powered colloids and the suspension 
was studied with optical microscopy. The measurements 
show that the active colloids are hotter than in the pas- 
sive limit with an effective temperature that increases as 
the square of the parameter that controls activation, the 
Peclet number, a dependence that is highly reminiscent 
of the /^-dependence we found with our simulations. 

This summary proves that there is an increasing ac- 
tivity in the study of FDT violations and their interpre- 
tation in terms of effective temperatures in biophysics 
problems. A comprehensive experimental study of the 
thermodynamic properties, as we tried to present in this 
work within one system is, however, still lacking. In our 
opinion the active colloidal suspension used in [4]| ap- 
pears as the most suitable experimental system where to 
deal with this kind of questions. 



B. Discussion 

Let us now critically review our findings. We presented 
a study of active matter made of semi-flexible filaments 
in interaction. In the passive limit the conservative inter- 
action parameters put the system in a liquid phase. The 
equilibrium polymer melt has, though, properties that 
are different from the ones found in the limits of com- 
pletely rigid or totally flexible polymers. For example, 
the individual chains are not Gaussian. 

Under the action of mild motor forces, the structure 
and dynamics are signiflcantly modifled. The analysis of 
the chain radius of gyration suggests that the polymers 
become more compact for increasing forces. This is con- 
sistent with the behavior of thermodynamic quantities 
that suggest crowding under increasing forces. On the 
other hand, the study of the static structure factors sug- 
gests that the melt gets more disordered -as in a passive 
sample in contact with hotter baths- and that chains be- 
come more folded when / increases. As a consequence, 
the dynamics gets faster under stronger forces. This is 
shown by the fact that the relaxation time, as extracted 
from the evolution of the incoherent scattering function, 
decreases, while the diffusion coefficient increases, for in- 
creasing /. 
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Next, we developed an in-depth study of the effective 
temperature of the active system, by using a number of 
independent measurements. The resuhs of ah these mea- 
surements in the semi-flexible polymer sample, as well as 
the outcome of similar measurements in the self-propelled 
point-particle model we analyzed in [43] , are summarized 
in Fig. [201 

The dependence of T^ff on / conforms to our intu- 
ition. We dealt with motors that act randomly - and 
ignore their effect on the energy landscape of the system. 
They have the effect of disordering the sample and, con- 
sequently, the stronger / the higher Teff. In their study 
of point-particle active matter Q, Shen and Wolynes 
argued that motors that act selectively, only when they 
contribute to decreasing the system's free-energy density, 
lead to out of equilibrium dynamics with lower effective 
temperature than the environmental one. The interpre- 
tation is that this kind of motor orders the sample more 
than the passive system would be in equilibrium at the 
working bath temperature. A similar inversion has been 
observed in relaxing passive systems driven from an or- 
dered initial equilibrium condition to a point in param- 
eter space in which the evolution slowly tends to reach 
a less ordered final equilibrium state [73|. In this case, 
the dynamic configurations are also more ordered than 
the target equilibrium ones and, as a consequence, their 
effective temperature is lower than the one of the bath. 

The results discussed in this paper and in Ref. [i^ give 
a different perspective on the dynamics of active matter 
than the coarse-grained hydrodynamic descriptions [Tsl - 
[TtI . While in the mode-coupling study in [la| a strong 
frequency dependence of the effective temperature was 
found, in our simulations the behavior observed is, in 
contrast, very similar to the one advocated in [2^ to be 
the one of interacting glassy-like systems evolving in a 
small entropy production limit. Indeed, the global behav- 



ior of Teff in the active matter polymer system resem- 
bles the one computed in a variety of atomic models in 
the glassy regime [11], sheared super-cooled liquids [69| . 
driven vortex systems [70{ or even vibrated granular mat- 
ter [Tlj . Independent measurements of Tef / yield consis- 
tent results in the atomic and molecular cases, as shown 
in Fig. [201 giving support to the idea that T^ff can be 
considered to be a thermodynamic object, at least in the 
weakly driven cases studied here (/ < 1). It would be in- 
teresting to test whether a similar phenomenology arises 
in lattice models of Vicsek-type [1, [l^ with various inde- 
pendent measurements of T^ff yielding the same result. 
On the experimental front, it would be welcome to ex- 
tend the analysis of Te// in the dilute active suspension 
of artificial swimmers studied in |41| or search for T^f f in 
artificial samples such as the ones recently designed for 
their tunability [82| . 

Finally, it would be interesting to explore the possibil- 
ities offered by the approach described in this work to 
clarify the extremely complex phenomenology of cell me- 
chanical stability t3|. In particular, one could try to give 
an answer to the question as to how mechanical proper- 
ties like elastic moduli change with motor activity and 
whether these changes can also be rationalized in terms 
of the effective temperature. 
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